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1. Introduction 

Recent advances in observational astronomy have allowed maps of the cosmic microwave 
background to be constructed in more detail than ever before [1]. The availability of 
such maps, with good noise properties and controlled foreground subtraction, offers 
the exciting prospect of studying the earliest ages in the evolution of our universe in a 
relatively direct way [21 [3]. 

In order to connect any theory of early universe physics with the cosmic microwave 
sky, one must make predictions for properties of the curvature perturbation, £ [3], which 
sets the initial condition for those fluctuations in the matter and radiation densities 
which we can now observe as temperature anisotropies. To carry out this programme 
effectively, one needs both a model for the relevant physics and an efficient calculational 
tool with which to make predictions. 

One of the most attractive models which has been proposed to describe the evolution 
of the early universe is inflation (for a review, see Ref. [5]). In an inflationary scenario 
the universe is supposed to have undergone a phase of accelerated expansion in the very 
distant past. During this phase each light scalar field acquires a fluctuation, which 
is close to scale invariance when the rate at which the universe expands is almost 
constant. These fields each contribute a proportion of the total energy density of the 
universe, and the relative importance of these contributions is sufficient to determine 
how ( is composed of the separate fluctuations in each light scalar [6j [7J [H [9j [10] . The 
fluctuations themselves can be computed, for example, by using the methods of quantum 
field theory. Inflation thus provides a framework within which the properties of £ can 
readily be calculated, and for a large class of models these predictions are in extremely 
good agreement with the observational data [llj. 

Which properties of ( should we compute? Any measurement of the CMB 
anisotropy can be framed in terms of the n-point expectation values of the curvature 
perturbation, and current data are primarily sensitive to the simplest such expectation 
value, which is given by the two-point function or power spectrum. Although this 
observable has sufficed to obtain a good deal of precise information about the very 
early universe, there are limits to how much we can learn from it. Therefore, in 
order to discriminate effectively between the different models of inflation, it will be 
necessary to broaden our observational methods so that higher n-point functions of £ 
become experimentally accessible, and at the same time develop theoretical predictions 
for these functions in all relevant models. On the experimental side this process has 
been underway since the Cosmic Background Observer (COBE) satellite made the first 
all-sky map of the CMB in the mid 1990s p2} [13]. Such efforts have borne fruit with the 
availability of more sensitive measurements taken by the NASA Wilkinson Microwave 
Anisotropy Probe with high angular resolution [H[[ll], and will be improved still further 
by ESA's Planck satellite, due for launch in mid-2008. The first year WMAP data release 
placed significant constraints on the level of non-gaussianity which could be present in 
the CMB (often quantified in terms of a so-called "non-linearity parameter," /nl [15j). 
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and this quantitative estimate has subsequently been refined [HHELTl, leading to a recent 
high-confidence exclusion of /nl = [18] . In the long-term, it may be possible to obtain 
even better estimates from observations, perhaps based on maps of the 21-cm emission 
of neutral hydrogen [191 [20]. O n the theoretical side, predictions for the non-gaussianity 
generated in a large collection of relevant models have become available over the last 
several years ^TJ 1221 (EH 123 1231 123 12SI 123 12SI 123 EDI IBl ED E21 E3J , following Maldacena's 
successful calculation of the bispectrum produced in single- field, slow-roll inflation [34]. 

These developments mean that non-gaussianity is now a standard cosmological 
observable, comparable to the spectral index n of the scalar power spectrum, and is a 
powerful discriminant between competing models. It is therefore extremely desirable to 
have at our disposal a simple calculational scheme — analogous to the familiar formula 
which allows the scalar spectral index to be estimated [5] — which can be used to predict 
the non-gaussianity generated in any model of our choice. 

Most computations of the primordial non-gaussianity have been carried out using 
the tools of quantum field theory, and have been explicitly framed in the context of 
the Lagrangian formalism, which leads to a method of computation closely related 
to many calculations in particle physics [341 l29[ 135] . For example, when expanded 
perturbatively, the Lagrangian formalism naturally gives rise to a variant of the Feynman 
diagrams which are familiar from the calculation of scattering amplitudes. There is 
no doubt that this is a useful development, which immediately permits methods and 
intuition developed in the context of quantum field theory to be imported into cosmology. 
Nevertheless, there is considerable merit in exploring alternative calculational strategies, 
either to enlarge the class of theories in which the non-gaussianity can be computed, or 
to take advantage of technical simplifications in the computation. 

In the present paper we return to the more familiar method of employing the field 
equations, which does not invoke a Lagrangian. This is how the tree-level spectrum 
(H/2ti) 2 of a scalar field in the de Sitter spacetime was first calculated [36], and also 
how the tree-level spectrum of the curvature perturbation was first calculated using the 
Mukhanov-Sasaki field equation [37J [38] More recently it has been shown how the 
method allows one to calculate all expectation values of field perturbation to all orders 
in de Sitter space [40] . Here we use the method to calculate the tree- level bispectrum of 
the inflaton field perturbation, to second order in the field perturbation and including 
the associated metric perturbation. 

The method of field equations permits us to make contact with the cosmological 
literature, which has traditionally approached perturbation theory from the standpoint 
of the Einstein equations. A second advantage is that it allows an immediate calculation 
of the non-gaussianity generated in any model which can be cast in the form of an 
effective Klein-Gordon equation, whether or not the model is descended from an effective 
action principle. We do not invoke the Einstein-Hilbert action, only the Einstein field 
equation together with the inflaton field equation which is obtained from the flat-space 

X "Tree-level" here means leading order in the perturbation, which indeed corresponds to the tree-level 
Feynman graph [391140]. 
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Klein-Gordon equation by the replacement rj^ — > (which is the content of the 
equivalence principle). We are able to give a simple, essentially "mechanical" process 
for computing the bispectrum in such a model. 

The outline of this paper is as follows. In £J2] we construct an explicit solution for 
the Heisenberg picture field directly from the equations of motion. In £J3] we use this 
solution to construct the three-point correlator of the field fluctuations. This is the 
essential ingredient in a calculation of the bispectum of £, since it serves as the required 
initial condition for the well-known 5N formula. In §H]we write down an explicit formula 
for the non-linearity parameter /nl before giving a brief summary of the calculation in 
§51 which concludes with a discussion. 

Throughout this paper, we work in natural units where the Planck mass is set equal 
to unity, Mp = (87rG)~ 1//2 = 1. The background metric is taken to be of Robertson- 
Walker form, 

ds 2 = -dt 2 + aitfdij dx i dx\ (1) 

where the scale factor a(t) obeys the usual Friedmann constraint 3H 2 = p, with p being 
the energy density. Throughout, H = a/ a is the Hubble parameter and an over dot 
denotes a derivative with respect to t. It is sometimes more convenient to work in terms 
of a conformal time variable, 77, which is defined by 77 = f°° dt'/a(t'). When 77 is used, 
it is useful to define the conformal Hubble rate 7i, which satisfies 7i = a' /a = aH. 

2. Constructing the Heisenberg-picture field 

For simplicity we will assume there is only a single scalar degree of freedom, although 
our method generalizes easily to the case of many scalar fields. The perturbed metric 
can be written 

ds 2 = -N 2 dt 2 + h i3 (dx l + N* dt) (dx j + N j dt) . (2) 

We work in the uniform curvature gauge which corresponds to setting hij equal to its 
unperturbed value [given by Eq. (CTJ], where we explicitly assume that all tensor modes 
are absent J§| The quantities N and N l are, respectively, the lapse function and the shift 
vector [4"T] . These are unambiguously determined by the Einstein constraint equations 
once a gauge has been chosen. 

Let be a canonically normalized field, whose background value sources the 
metric (CQ). We assume that <fi is driven by a potential V, which supports inflation 
in some region of field space. One then allows to be perturbed by a small amount 5(f), 
which obeys an effective Klein-Gordon equation. Perturbation theory is an expansion 
in powers of 5(f). 

In the cosmological literature a further decomposition is sometimes made (see, for 
example, Ref. [12]), by separating some part of 5(f) which can be called 5(f>\. One then 

§ A background of gravitational waves would not contribute to the bispectrum to leading order in the 
perturbations; its contribution would therefore be strictly subleading. For this reason, the truncation 
to zero tensor modes should be acceptable for the purpose of computing /nl- 
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writes all quantities in the perturbation theory, including 50, as an expansion in powers 
of 501, 

50 = 50x + ^50 2 + ■ ■ • + —5<p n + ■■■. (3) 

It is necessary to impose some arbitrary auxiliary condition to determine 50i, which 
is usually done implicitly by demanding that 50i obey a linear equation of motion. It 
follows that 50i exhibits precisely gaussian statistics. This formulation of perturbation 
theory is most useful when 50 becomes increasingly gaussian at early times, which 
implies that the 50 n for n ^ 2 vanish at past infinity. These two pictures of perturbation 
theory have complementary merits, and can be used interchangeably. 

We will work in the second picture and choose 50i to be that part of 50 which obeys 
precisely gaussian statistics, as discussed above, and formulate perturbation theory as 
an expansion in powers of 50i. As we have described this implies that 50i should be 
chosen to obey a linear equation of motion, so that each Fourier mode evolves according 
to the equation 

50" + 2^50; + £; 2 50i = 0, (4) 

where a prime ' denotes a derivative with respect to rj. Eventually we are aiming to 
calculate only to leading order in the slow-roll expansion and therefore higher-order slow- 
roll contributions in Eq. (j3]) have been omitted. It is useful to introduce the notation 
9k for solutions of (J3J, so that 50i(k, rf) = 6k{i])- 

So far, this theory is entirely classical. Quantization of 50i is straightforward, and 
proceeds by the canonical method (see, eg., Ref. [S]). The quantum field corresponding 
to a conjugate pair of solutions {6k, 6k} is 

50i(x,r / ) = J -0- 3 ^{at6 k (v) + a-iAfo)} , (5) 

where the normalization of the {6k, 6k} h as been adjusted so that the canonical 
commutation relation for 50 and its conjugate momentum is obtained when {a^, a k } 
obey the usual creation-annihilation algebra, 

[a k ,a*,] = (27r) 3 5(k-k'). (6) 

On the other hand, 502 obeys a non-linear equation with quadratic source terms. 
In this paper we are going to use the slow-roll approximation to control these sources. 
As usual we take slow-roll to be defined by the following set of conditions: 

H 1 fV'\ 2 1 $ , 
V" 

r] = — ~ — -^ + e<l. 8 

V H<j> 

(In particular, one should take care that rj is not confused with the conformal time.) 
The non-linear equation which describes 502 can now be written in the form 

50 2 ' + 2H502 + k 2 5<p 2 = -a 2 V"'S(f>l + F 2 (50i) + G 2 (50'i), (9) 
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where V" is the third derivative of the potential for 0, and F 2 and G 2 are source terms, 
which are respectively quadratic in 8<pi and 5(j)' x . It has recently been shown that [HIH5] 



where we have dropped terms that are suppressed by factors of e and 77. It is useful to 
separate the a 2 V" term from the remainder, because its different dependence on the 
scale factor endows it with a distinctive time dependence which we shall see below leads 
to a logarithmic divergence at late times. 

Among the small terms — subleading in powers of e and 77 — which have been dropped 
in Eqs. (l9l)- (fTT]) are source terms proportional to the cross-product 501(50^. Such terms 
also lead to logarithmic divergences which have the rough form ~ ej N, where for any 
mode of wavelength k the quantity iV = In |fc?7| measures by how many e- foldings it is 
outside the horizon at some time of observation, 77, and a subscript V denotes evaluation 
at horizon crossing. Since we are assuming that slow-roll applies at that time these terms 
are small when N w 0, but they grow in magnitude with N and when iV ~ 1/e* they 
can no longer be ignored. Indeed, such terms form part of a tower of divergences which 
have the rough form ~ eT +1 ^ 2 N m for m ^ 1. It follows that the slow-roll approximation 
will be close to breaking down if we wish to evaluate expectation values at the end of 
inflation, where the scales of interest are N e-foldings outside the horizon (with values 
of order N ~ 60 being typical), but in many theories e* ~ 1/N. This would lead to a 
nonsensical perturbation theory in powers of unity. 

This argument shows that we cannot use an expansion in powers of slow-roll 
parameters to predict expectation values at the end of inflation. One should instead 
think of these divergences as terms in a Taylor series expansion around the time of 
horizon exit, where N = [46] . In principle they could be resummed by the method 
of the renormalization group or an equivalent technique, after which they would merely 
correspond to the classical time evolution [171 HH1 SSI ED]. However, there is no need 
to perform such a complicated resummation. Instead, we believe it is most accurate to 
evaluate all expectation values almost immediately after horizon exit, where slow-roll is 
an excellent approximation and all strictly positive powers of N are entirely negligible. 
One must then use some other method, which does not rely on an expansion in terms of 
e-foldings since horizon exit, to evolve these expectation values until the desired time of 
evaluation. This is equivalent to the argument of the renormalization group, but makes 
use of the known simplification that evolution outside the horizon is simply classical. 
For the curvature perturbation the 5N formula is very convenient 0, [361 13 113 HO- 
More generally one can use the separate universe approach [511 [7J [H [251 E2j, or an 
equivalent gradient expansion [HUE]- It may also be possible to use some formulations 
of conventional perturbation theory, provided they avoid the appearance of powers of 
N [S] . The crucial point is that if the correlation functions are evaluated immediately 



and 



F 2 = - 



G 2 




(10) 



(11) 
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after horizon exit, it is only necessary to compute the leading term in the Taylor series, 
which will be the constant term in any model giving rise to a Klein-Gordon equation 
of the form ([9]). This term is given by the lowest-order slow- roll approximation. 

In theories which are more general than Einstein gravity coupled to a scalar field, 
one might encounter examples where the constant term in the Taylor expansion is absent, 
or where the linear term has a comparable magnitude. An example of such a model 
is the case of non-local inflation studied by Barnaby & Cline [53] , and we shall see 
below that it can occur even in single-field, slow-roll inflation if the third derivative 
of the potential is exceptionally large in comparison with the first and second. In 
such cases one can compute enough terms in the Taylor series to obtain a satisfactory 
approximation before using the separate universe picture or some equivalent technique 
to obtain the superhorizon evolution. If any 5</>i5</4 terms participate in this process, 
then they can be computed using a small modification of the argument outlined below. 

In the classical theory, Eq. ([9]) can be solved by the use of a Green's function. The 
correct choice is the so-called causal or retarded Green's function, which in momentum 
space satisfies 

{0 n <r t 

where {6 k , 9 k } are a complex conjugate pair of solutions to the non-interacting Klein- 
Gordon equation (j3J) for comoving wavenumber k, normalized as in Eqs. ([5])-(j6]), and a 
bar denotes complex conjugation. To use Gr^ to solve for 6(p 2 , one must translate the 
source terms F 2 and G 2 into Fourier modes. The behaviour of each interacting mode 
can then be computed, and the result reassembled, giving 

6Mv > x) = / (0 eiqx { /I dr / Gr<?(r? ' r)5(q " kl " k2)l52 } ' (13) 

where S 2 is defined by 

S 2 ee [-a\r)V" + T 2 {\^M)} ^i(k 1 ,r)^(k 3 ,r)+g a (k 1 ,k a ) d<y ^ 1 ' T) <M ^' t) > (14) 
and the momentum space source functions T 2 and Q 2 satisfy 

- ^ ■ k 2 + 2k\ + 1 [kl4 + foX ■ k,]) (15) 

fes -Ji^ + k 1 ' k2 -WTW [kl + kl - k2] )- (16) 

In the quantum theory the same construction applies, with the understanding 
that we are to substitute the quantum field 6<f>i(q) ee a^9 q (r]) + a_ q ^ g (r]) for the non- 
interacting solution 5(f>i (q) = 6 q in the source term S 2 . Thus, once the quantization of 
54>i has been determined, this is sufficient to fix the quantization of the S(j) n for n ^ 2. 
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3. The three-point correlator 

The final step is to use our solution for the interacting Heisenberg field to compute 
the three-point expectation value, (50(k 1 )50(k 2 )50(k 3 )). Since 5(fii obeys precisely 
gaussian statistics, it is forbidden to have a non-trivial three-point correlator. The 
leading contribution therefore comes from a correlation between the second-order part 
of one of the fields with the first-order part of the remaining two, so that the correlation 
is schematically of the form (<5</>i<50i|<50 2 ) ~ h(S<piS<pi8<pi * 6~(fii), where V denotes a 
convolution. The resulting four-<50i expectation can be evaluated by use of Wick's 
theorem, together with the Wightman function for S(j)i, 

(5<Mk,7/)5<Mk',r/)> = (27r) 3 5(k + k%(r/)^(r/). (17) 

Notice that there is no time ordering in Eq. (fTTj) . 



3.1. The V" term 

Consider first the V" term, leading to a logarithmic divergence as described above. The 
three-point expectation value sourced by a term of this type was first calculated by 
Falk, Rangarajan and Srednicki [M], and was later reconsidered by Zaldarriaga, who 
calculated it in the context of the inhomogeneous decay rate model [47J. 

This term gives a contribution to the three-point function which can be written 

(50(k 1 )50(k 2 )<50(k 3 )) D - i(2vr) 3 5(k 1 + k 2 + k 3 ) f dr a(r) 4 V w x 

J — oo 

K(r% 3 (rj) - 9 k3 (r)0kM] M^M^M^W + 

[0Ur)hM - h x {r)e k M] e k MhMK{r)h 3 {r)\- (18) 



It is easy to verify that four of the six terms in (1181) cancel against each other, leaving 
the sum of a quantity and its complex conjugate, 

(<ty(ki)<ty(k 2 )<ty(k3)> D - i(27r) 3 5(k! + k 2 + k 3 ) ^ dTa(T) 4 V'"l[9 h (T)6 k M 

J -co i 

+ complex conjugate. (19) 

where i G {1,2,3}. Such cancellation always occurs in the Lagrangian version of the 
in-in formalism, leading to expectation values which are manifestly real. In the present 
formalism any cancellation is contingent on the possibility of factoring the 8, 9 terms. 
As we shall see below this does not always occur, leading to calculations which are 
somewhat longer than the Lagrangian equivalent and to the loss of manifest reality. At 
the end of the calculation, of course, the two methods must produce equivalent results. 

The integral in f[T9l is badly divergent, but almost all the divergent terms are 
imaginary and cancel out of the final result. This is easiest to see after integration by 
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Tr2ym 

(50(k 1 )^(k 2 )^(k 3 )} D (2 7 r) 3 5(k 1 + k 2 + k 3 )^^3 x 

-lkl + k t Y,hk 3 + \{\ + l + ^\hnA}Y, k U > ( 2 °) 

i<j i / 

where i G {1,2,3}, k t = J^^ki, and for the purposes of obtaining an analytic solution 
we have assumed that the ki have some approximately equal magnitude, which defines 
a more or less unique time of horizon exit 77*. A subscript V denotes evaluation at this 
time, and 7 ~ 0.577216 is Euler's constant. 



3.2. Zero-derivative terms 

Now let us focus on those quadratic source terms in Eq. (J9]) which enter with no time 
derivatives, described by the function F 2 . Such terms make a contribution to the three- 
point expectation value of the form 

(50(k 1 )50(k 2 )50(k 3 )) D i(27r) 3 (5(k 1 + k 2 + k 3 ) dr a(r) 2 x 

J —00 

|f 2 (a ;i , k 2 - h) [e k3 (r%M - e k3 {r)e k M] hMhMKij)B k2 {r) + 

F 2 (*i, h; k 2 ) [e k2 {r)d k M - h 2 {r)0kM] k MK(v)0 kl (r% 3 (r) + 

¥ 2 (k 2 ,k 3] k 1 ) [e kl {r)e k M-hAr)e k Av)\ d k MdkMh 2 (r)e k3 (r)^, (21) 

where the fields which participate in the expectation value are each evaluated at time 
r/*. The quantity F 2 (a, 6; c) is to be obtained by symmetrizing JF 2 (a, b) between a and b 
with weight unity. (This follows from the factor of 1/2 which is used in the definition of 
S(p 2 ; if this factor is absent, the the relative weighting of F 2 must be adjusted to match.) 
However, the result is subject to the constraint a + b + c = and therefore can be 
written in several equivalent forms. For this reason it is useful to fix the representation 
unambiguously, following Maldacena [34], by using the relation between {a, b, c} to 
eliminate all dot products between distinct vectors, leaving a result which depends only 
on the magnitudes {a, b, c}. We find 

UK h; h) = -I (^(kl + kl) - _ M) . (22) 

If F 2 (a, b;c) were symmetric over permutations of {a,b,c}, then as above the six 
terms in Eq. (12 ID would partially cancel among themselves. Owing to its construction, 
however, F 2 is in general symmetric only under interchange of a and b and fails to exhibit 
any symmetry under exchange of c with {a, b}. Therefore, it is not obvious that Eq. ( 12~TI) 
is real. To go further and demonstrate equality with the Lagrangian formulation, it is 
necessary to compute each of the six integrals. 

1 1 Note that this result differs from that given by Zaldarriaga by a sign and the replacement 7 > 7+ 1/3. 
We thank Filippo Vernizzi for pointing this out to us. 
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Once this is done, we find that the zero-derivative contribution to the expectation 
value becomes 

(50(k 1 )^(k 2 )^(k 3 )} D (2vr) 3 5(k 1 + k 2 + kg) x 

" ;/ ( fi¥ 2 (k u k 2 ; k 3 ) + f 2 ¥ 2 (k u k 3 ; k 2 ) + / 3 F 2 (A; 2 , h; (23) 



8 EL kf 

where the momentum factors fi obey 



2ki(kf + Ahk 2 + k 2 2 - kf, 

(h + k 2 - fc 3 ) 2 fc? 
2k\(k\ - Ak^ + kj- k\) 



{k\ - (hi - k 3/ 



2\2 



2kl(k% + 4k 2 k 3 + ifc| - Jfc^ 



(A* + £; 2 - A; 3 ) 2 A; 



and we have again assumed that the ki have an approximate common value, k. The 
expectation value is to be evaluated at a time just after horizon exit of the mode with 
wavenumber k. As above, H* denotes the value of the Hubble parameter at this time. 
Note that the fi do not vary from theory to theory; they are fixed by the structure of 
the interaction, and once calculated do not need to be calculated again. The different 
theories which may govern the interactions of 5(f) influence the three-point expectation 
only through the term ¥ 2 (a,b;c). Similar remarks apply for the n-point correlation 
functions for all n. 

One might have expected that the fi should be related by cyclic permutations 
among the various fcj. In fact there are extra signs which are introduced, because the fi 
arise from integrating over combinations of the mode functions 9 ki and their conjugates 
6 ki . This explains the apparent lack of symmetry in Eq. (J23J- 

3.3. Two- derivative terms 

We can carry out a similar calculation as in the previous section for the terms in the 
5(j) 2 equation which are quadratic in S<pi, given in (fTTj) . One arrives at a contribution to 
the three-point expectation which can be written 

(<50(k 1 )50(k 2 )50(k 3 )) D i(2vr) 3 5(k 1 + k 2 + k 3 ) / * dr a(r) 2 x 

J — oo 

G 2 (h, k 2 - k 3 ) [9 k3 (r%M - 6 k3 (r)e k3 (r])] hM^M ^^ + 

G 2 (h, h- k 2 ) [d k2 {r%M - h 2 {r)e k M] hMd k M d6k ^ T) d9k ^ T) + 
G 2 (k 2 ,h;h) [e kl (r)e k M - flfafrOMq)] 0kM0kM d6k : ir) d ^ (r) ), (25) 



where G 2 (a, b; c) is obtained from Q 2 (a, b) by a rule analogous to that relating F 2 and 
T 2 . one symmetrizes over a and b, and then uses the relation a + b + c = 0to eliminate 
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G2(*i, fc 2 ; jfe 3 ) = ttj o? 2?2 3 - (26) 



dot products of distinct vectors. This gives 

g (fc? + fc 2 2 )(fc; + A: 2 - A; 2 , 

In Eq. there is no cancellation between terms even in the case that 62(0, 6, c) 
is totally symmetric. However, we can follow a similar line of argument used in the case 
of the zero-derivative terms and find 

(^(k 1 )50(k 2 )50(k 3 )) D (27T) 3 5(k 1 + k 2 + k 3 ) x 

Hi 



8 EL kf 

The two-derivative momentum factors gi are defined by 

4 a*? 



giG 2 (h, k 2 ; k 3 ) + g 2 G 2 (h, £3; fe) + g 3 G 2 (k 2 , k 3 ; h) \. (27) 



+ k 2 - k z fk 2 t ' 

{k\ - (^ - fc 3 ) 



9i 

^ 3 = (kf - (k 2 + k 3 yy ' (28) 

We can now assemble the zero- and two-derivative terms to obtain an overall 
contribution to the three-point expectation value which corresponds to 

(50(k 1 )^(k 2 )50(k 3 )) D (2 7 r) 3 5(k 1 + k 2 + k 3 ) x 

in complete agreement with previous calculations using the action-based approach 
[23 [551. 



4. The non-linearity parameter /nl 

The three-point correlator of 50 is not itself the object of primary interest, since it 
is not directly observable in the cosmic microwave background. Instead, experiments 
conventionally set limits on a parameter /nl which is related to the three-point correlator 
of the curvature perturbation, (. 
We define /nl by 

(C(k 1 )C(k 2 )C(k 3 )) = (27r) 3 <5(k 1 + k 2 + k 3 )^ NL ^P ? (^)Pc(A;,). (30) 

i<j 

With this choice of sign and prefactor, the parameter /nl corresponds to the one 
conventionally used in the analysis of CMB observations [56], although other sign 
conventions exist |34j . 

For a general theory one can use the 5N formula to show that 

dN 1 d 2 N 

C(x) = SN = %t-6Mx) + 2 V^* (X) + " ' ' (31) 
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which holds in coordinate space and is easily generalized to the case of multiple fields. In 
this equation, the <50 are evaluated on a flat hypersurface at time 77*, and N is a function 
of both 77* and the time of observation. To second order in 50, ( can be calculated in 
terms of the so-called "separate universe approach," and takes the fornjjjl 

C(k) = y^^(k) + ^ ( X " I7) / 5 M±i - q)^(q) + • • • , (32) 

where denotes terms of higher order in 50* which have been omitted. In this 
expression, * denotes evaluation at some time t after the mode with wavenumber k has 
left the horizon, after which £(k) maintains a constant value. 

To complete the calculation of /nl we will use two more slow-roll parameters: 
V'V" 

(33) 



V 2 



a 



V 3 

and the slow-roll relations 



<yo 2 (d 4 y/d0 4 ) 

Latic 
d(ln e* 



(34) 



and 



diV* 

div; : 

dC 2 



~ 4e* - 2?7*, (35) 
2t*V*-£, (36) 



(1V -4e4i-V*C-<- (37) 

The slow-roll approximation, which we defined through Eqs. (CI) — (JHJ) , does not in 
itself place conditions on £ 2 and a 3 . Barring rapid oscillation of these quantities as 
functions of 0, or a narrow spike, we certainly need them to be much less than unity to 
preserve \r]\ <C 1 over many e-folds. At a generic point it is reasonable to expect [57] 
|cr 3 | -C |£ 2 | <C 1 77 1, but that will obviously fail if 7/ goes through zero. 
Using (|32|) to compute (CCC) an d using Eq. (I30I) for / NL gives 

g Jnl = -e* - V* + L I 3 + 7 + 

where we recall that iV* measures the number of e-folds which elapse between horizon 
exit of the mode k and the time of evaluation, 77*. 

If taken at face value, Eq. fl38|) may suggest that we can obtain an / NL as large as 
we please by allowing a large number of e-foldings outside the horizon, making N* very 
large. However, that cannot be the case because the separate universe approach shows 
that C is conserved after horizon exit during single-field inflation. (More generally, ( is 

% In writing this expression, we are assuming dV/d<j) > which can be accommodated by reversing 
the sign of 4> if necessary. 
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conserved outside the horizon if isocurvature perturbations are negligible so that there is 
a unique relation between pressure and energy density [581 [H [59].) From fl36|) we see that 
the explicit iV* dependence of /nl is cancelled by a contribution to the iV* dependence of 
77*. The remaining time dependence of /nl involves contributions proportional to e*?]*, 
rjl and of, which would be cancelled if we took the calculation to a higher order in the 
field perturbation and abandoned the slow-roll approximation. 

The present example is trivial because ( is time- independent. However, a similar 
caution applies to multiple-field inflation models where ( is time-dependent. There also, 
powers of iV* may also appear in the expression for /nl- If these terms are to source 
large non-gaussianity (which occurs, for example, in the curvaton scenario [601 1521 f3Tj ) . 
then they should be handled using a non-perturbative approach rather than relying on 
"naive" expressions such as Eq. ( !38l) . 

5. Conclusions 

In this paper we have shown that the bispectrum of an interacting but canonically 
normalized scalar field can be calculated during an inflationary epoch using the details 
of its field equation directly, without constructing an effective action and using the rules 
of the in-in formalism. 

Where a Lagrangian formulation exists, the action and the field equations to which 
it gives rise are manifestly equivalent. In this point of principle, our calculation 

is a straightforward rewriting of the usual one, although it may prove to be more 
convenient in examples where one works with the field equation from the outset. On 
the other hand, in models where no Langrangian formuation exists, or none is known, 
our formula enables the non-gaussianity to be computed easily and compared with 
experiment. 

The recipe is straightforward. For reference, we summarise the steps here. 

• Write the model in terms of a scalar field 0, and determine the field equation for 
the perturbation in this field, 5(f), defined in the uniform curvature gauge. 

• Separate d~<j) h^o a term 5(f>i which obeys a linear equation of motion, and a part 
S4>2 which is quadratic in 8<pi. Write down the evolution equation for 50 2 . 

• Determine the Fourier coefficients JF 2 and Q2 (respectively) of the terms quadratic 
in S(f>i and S<p' v 

• Construct the quantities F 2 and G 2 from JF 2 and Q 2 - 

• Insert F 2 and G 2 in Eqs. ( |23l) and ( 1271) . respectively. Simplify the resulting 
expressions to give the complete three-point expectation value. 

We have stressed that the computation of the 5(f) expectation value soon after 
horizon crossing is a simple exercise. However, we wish to emphasize that by itself this is 
insufficient to obtain a prediction for the non-gaussianity which is observed in the CMB. 
For that, one needs to make a prediction for the expectation values of the curvature 
perturbation, (, using the the classical evolution on superhorizon scales and the 5(f) 
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expectation values evaluated soon after horizon crossing as an initial condition. This 
process avoids the appearance of large logarithms which would require resummation. 

In the context of Einstein gravity coupled to a single scalar field this is simple 
to implement. The separate universe picture can be used in conjunction with 
rigorous results concerning the conservation of £ on superhorizon scales to make robust 
predictions, especially in the case where isocurvature modes are absent. In more general 
theories it may be necessary to augment this procedure or find a replacement in order 
to make reliable estimates of the non-gaussianity which can actually be observed in the 
CMB. 
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